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ABSTRACT 


A number of oscillatory phenomena occurring in 
lectrical circuits are treated in a series of three 
ulletins. In this first bulletin an introduction to 
he theory of almost sinusoidal oscillations in non- 
inear systems, its problems and methods, is given 
ith particular emphasis on simultaneous oscillations. 
he other two bulletins treat the problems of syn- 
hronization and transient phenomena. 

All of these problems are handled by a simple uni- 
orm method in which the nonlinear element, such as a 
acuum tube, is replaced by an "equivalent" linear 
mpedance. The value of this linear impedance is 
hosen so that, in the first approximation, the be- 
avior of the oscillatory system will remain the same. 
nce the circuit has been linearized methods well 
nown to the electrical engineer can be applied to it.’ 

This procedure is a direct expansion of the theory 
f "equivalent linearization" developed by Kryloff and 
ogoliuboff. ‘2)3)* The results obtained by it are 
dentical with those calculated by more complicated 
ethods. 


superscript, parenthesized numbers refer to corres-— 
sondingly numbered entries in the bibliography. 
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1. INTRODUCTION 


A great many of the phenomena occurring in nature 
which lead to nonlinear differential equations can be 
adequately described by linear differential equations. 
To arrive at such a linear differential equation, it 
1s necessary to neglect in the physical system those 
factors which have negligible influence on its be- 
havior. The desirability for linearization arises 
from the fact that the solution of linear differential 
equations is relatively easy since extensive mathemat- 
ical apparatus is available for their treatment; for 
nonlinear differential equations, on the other hand, 
1t 1s almost completely missing. The validity of such 
a linearization depends, of course, on the exact be- 
havior of the phenomenon under consideration. Thus 
for a triode oscillator the frequency can be found 
readily and with good accuracy from the linear har- 
monic differential .equation which is obtained by neg- 
lecting the nonlinear characteristic of the triode; 
but this nonlinearity cannot be neglected if the am- 
plitude of oscillation is being found. 

There is in nature a considerable variety of pheno- 
mena that cannot be expressed adequately by linear 
differential equations; that is, some important aspect 
of the behavior would be lost through such a lineari- 
zation. Some of these phenomena are, for example, 
turbulence in liquids, the generation of subharmonics 
by an inductance with an iron core, and the virtual 
cathode in a diode with free space-charge. 

Of this variety of phenomena leading to nonlinear 
differential equations those concerning oscillations 
are of great importance in many fields of science and 
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engineering, and their study is of considerable signi- 
ficance. 

Nonlinear mechanics deals almost exclusively with 
oscillatory phenomena. The reader who is not familiar 
with the subject is referred to such excellent books 
as Minorsky’s Introduction to Non-Linear Mechanics, 
Andronow and Chaikin’s Theory of Oscillations, and 
Stoker’s Nonlinear Vibrations. The field of nonlinear 
mechanics today is still fairly new. It has to op- 
erate under a severe handicap, inasmuch as an exact 
mathematical background is missing. Because of this 
lack of rigorous mathematical justification, the vali- 
dity of some of the methods is questionable. 

An important distinction in nonlinear mechanics is 
that made between conservative and nonconservative 
systems. Conservative systems contain a nonlinear 
storage element - for example, an inductance with an 
iron core. Nonconservative nonlinear systems, on the 
other hand, contain a nonlinear damping element, such 
as an electron tube; this bulletin deals exclusively 
with problems of the second type. 
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Figure 1.1 The Methods of Nonlinear Mechanics 
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|. The Methods of Nonlinear Mechanics 


No single general mathematical method is available 
which permits a treatment of all oscillatory problems. 
The cases that can be treated by the different methods 
are shown in Fig. 1.1. The topological method is 
best suited for the treatment of autonomous systems 
having any degree of nonlinearity but only one degree 
of freedom. The analytical methods can be applied to 
systems having a great amount of nonlinearity and one 
degree of freedom or to nearly linear systems with one 
or several degrees of freedom. A somewhat wider range 
of problems can be solved by methods involving analog 
or digital computers, but up to the present time no 
great use has been made of this possibility. 

In the topological methods the trajectories corres- 
ponding to the solutions of the nonlinear differential 
equations are studied in a two-dimensional system, the 
"phase plane," with the dependent variable as the ab- 
scissa and its time-derivative as the ordinate. °*7") 
For example, the solutions of the nonlinear differen- 
tial equation 


y) 
Soe i ene ee = 80 
dt2 dt (1.1) 


can be represented as trajectories in a plane withu 
and du/dt=w as the coordinates. Equation(1.1) 1s 
identical with 


dw. -u2) +U= 
u(1-u“) - 0. (1.2) 


[fhe solutions of this equation can be found graphi- 
sally by the method of isoclines. The solid-line 
“urves of Figs. 12.1 and 12.2 correspond to the peri- 
»ndic solutions of Equation (1.1) for two different 
values of the parameter Lt. The topological methods 
result, in some cases, in definite statements about 
che existence and stability of periodic solutions. 

The three analytical methods in use at present are 
named after their originators - Poincare, van der Pol, 
and Kryloff and Bogoliuboff. Not essentially differ- 
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ent, they give identical first approximations to the 
solutions of nonlinear differential equations. For 
any of the three methods very little is known about 
the validity of the higher approximations, and reason- 
able doubt as to their convergence is justified. 


2. The Method of Equivalent Linearization 


The method chosen for the discussion in this bul- 
letin is an analytic one developed by Kryloff and 
Bogoliuboff called the method of "equivalent linear- 
ization.” '43) Since this method has been developed 
primarily for the use of electrical engineers, the 
scope of the investigations reported in this bulletin 
is limited to the electrical circuit, and the method 
1s presented in the form of a "cook book recipe." 
Rather than to attempt any mathematical justification, 
the validity of the method is verified by agreement 
between the results obtained by it and experimental 
results. 

From physical considerations it can usually be 
assumed that in a given nonlinear circuit with high 
Q’s, either the currents or the voltages are nearly 
Sinusoidal, as for example the voltage v in Fig. 
3.2. By the method of equivalent linearization the 
almost sinusoidal current or voltage is replaced by 
a purely sinusoidal one. If the content of the 
higher harmonics in the original is small, then this 
substitution does not lead to appreciable errors. 

If it is assumed that some of the variables are 
purely sinusoidal, then the other contains the higher 
harmonics. (For example the current i in Fig. 3.2.) 
Since it can be shown that the behavior of the non-— 
linear circuit is defined by the real and imaginary 
average power generated by the nonlinear element and 
delivered to the linear elements and since this power 
is independent of the higher harmonics, they can again 
be neglected. Because under these assumptions the 
current flowing through the nonlinear element and the 
voltage across it are both sinusoidal, it can be re- 
placed by an "equivalent" linear impedance in such a 
way that the behavior of the system remains essen- 
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tially unchanged. Once this equivalent linear im- 
oedance is found, the theory of linear networks, well 
known to the electrical engineer, can be applied to 
the system. Thus the main advantage of this method is 
that it makes use of familiar concepts. To the mathe- 
natician, not trained to think in terms of physical 
circuits, it offers no special advantage. The method 
1s also well suited for the classroom and might well 
be incorporated into an advanced course in communica- 
tion engineering. 


Figure 3.1 The Standard Nonlinear Network 


3. The Nonlinear Circuits 


Since it is not practicable to treat all the pos- 
sible nonlinear oscillatory circuits separately, a 
standard network has been selected and is shown in 
‘ig. 3.1. This circuit shows a number of oscillatory 
>ircuits and constant voltage generators producing a 
sinusoidal voltage in series with the nonlinear ele- 
vent. If all the oscillators and all the oscillatory 


Figure 3.2 Oscillatory Circuit with One Degree of Freedom 
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circuits except one are omitted, there remains a cir- 
cuit with one degree of freedom corresponding to an 
ordinary triode oscillator as shown in Fig. 3.2. By 
omitting less elements, other cases arise. 


4. Conclusions 


Although the circuit shown in Fig. 3.2 is parti- 
cularly well suited for the theory of equivalent lin- 
earization, it severely restricts the number of prac- 
tical oscillators to which the following considera- 
tions can be applied directly; it excludes, for 
example, all networks containing several nonlinear 
elements. The method of equivalent linearization is 
not limited, however, to the circuit shown in Fig. 3.2 
but can be equally well applied to any circuit that 
has a small degree of nonlinearity. 


ll. SERIES EXPRESSIONS FOR NONLINEAR ELEMENTS 


General considerations of nonlinear elements are 
oresented, pertinent to the theory of equivalent lin- 
2arization developed by Kryloff and Bogoliuboff. If 
she nonlinear element can be approximated by a finite 
series of analytic functions, then the "linearized im- 
9edances" determining the behavior of the circuit can 
oe calculated rapidly by an operational method. 


5. Introduction 


In the literature on nonlinear mechanics, little 
consideration has been given to the nonlinear element 
itself. It is usually defined as "an element that is 
10ot linear." This definition is so vague and covers 
so many different elements that no general method of 
attack is possible. In this chapter a subgroup is 
selected to which the theory of equivalent lineariza- 
sion can be applied. For the other methods used in 
che theory of oscillators, for example, for those 
ising the phase plane, other definitions and subgroups 
ire more appropriate. It should also be pointed out 
shat the definition used here differs from the more 
sonventional one inasmuch as it describes only the 
-wo-terminal behavior of the nonlinear element and not 
_ts internal mechanism. Because linear elements are a 
legenerate case of nonlinear elements, everything said 
vere about nonlinear elements holds for linear ele- 
rents also. 

In the following discussion, consideration is re- 
stricted to nonlinear elements having the property 
hat a force or a voltage with a finite number of fre- 
juency components 


i 
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n 
= ‘oan Ae. J  ttO- 
W 2 zicos(wit 03) 


she | (Dees 
applied to two points of the element produces unique- 
ly a displacement or current which can be expressed by 
a uniformly almost periodic tunction of time and has 


components at the frequencies W lor our purpose we 


lp 
need only know that the uniformly convergent series 


Le 9) 


u(t) = 2 ajcos(Witto;) 


her (523 


falls anto the class of uniformly almost periodic 
functions, ‘2! 

This subgroup contains all elements usually en- 
countered in almost sinusoidal oscillatory problems. 
There are, however, a few notable exceptions where the 
"history" of the element has to be considered; for ex- 
ample, iron core or temperature dependent impedances, 

If the nonlinear element has time-varying para- 
meters, it 18 necessary to include corresponding terms 
in the right-hand side of Equation (5.1). If the in- 
side mechanism of the element is known, then the dis- 
placements or currents can be calculated when the 
force or voltage is given; otherwise, they have to be- 
measured, 

Some general considerations regarding nonlinear 
elements falling into this subgroup are presented, 
These considerations apply only to the case in which 
the nonlinear element forms part of an oscillatory 
eireuit and in which the energy stored or changed into 
other forms inside the nonlinear element is small re- 
lative to that stored in the linear parts of the cir- 
cult, In this case the oscillations are almost sinus- 
oidal, 


6. General Remarks 


In this chapter only nonlinear elements of th 
voltage-controlled type are treated, but all the con 
siderations discussed here also apply to current 
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-ontrolled elements and to mechanical or acoustical 
-lements as well. 


Figure 6.1. Example of an Oscillatory Circuit 


As an example Fig. 6.1 shows an oscillatory cir- 
uit composed of two resonant circuits, an external 
oltage E, and a nonlinear element NL. The two reson- 
nt circuits are formed by the inductances Lp, and Loo 


nd the capacitances Cy; and Cy9. The resonant fre- 


uencies of these circuits are Wy, and WO. 


Lt Oe 
Wo] ~ ae and Woo * 
V 191 Loong. 


t 1s assumed that the external voltage E is sinus- 
idal and has a frequency w3. The instantaneous value 


f the voltage across the nonlinear element is v, and 
he current passing through it is i. It is assumed 
hat the maximum value of the current i is small com- 
jared to the amplitudes of current in the resonant 
reuits. 

In the circuit of Fig. 6.1, the voltage can be 
oproximated by a finite trigonometric sum 


n 
wete e V; cos( ; tto, ). (6.1) 
ial! 


ne angular frequencies Wj are referred to in the 


ollowing discussion as the fundamental frequencies. 
ney are approximately equal to the resonant frequen- 


It is assumed that the Wj are not close to- 
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A voltage-controlled element falls into the sub- 
group discussed in Section 5 if the current can be 
expressed by a uniformly convergent trigonometric 
series, 


n 
1 * Io + 2 Ip jcos(witto;) 1 I, jsin(w;tto;) 
S=Veguarn re 
+ currents at combination frequencies 
and higher harmonics, (6.2) 


Some combination frequency or higher harmonics may be 
equal to one of the frequencies W;, and is then auto- 
matically aneluded in the second or third term of 
Nquation (6,2). 


Figure 7,4. The Equivalent Linear Network 


7. The Equivalent Impedances 


According to the theory of equivalent lineariza- 
tion the nonlinear element can be replaced by a linear 
element that produces a current 


n 
1 tea) + I,jsin(wjtto;). (7.19 


Such an element, a conductance G in parallel with ¢ 
capacitance C, is shown in Fig. 7.1. The values fo 
G and CG depend on the voltage v and are different fo 
each frequency. From Equations (6,2) and (7.1) thi 


values of G and C for frequencies w; are 
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ies I, 
eV epee 
i 
oo apm 
BEY (7.3) 


i; and C; are the "equivalent" conductance and capaci- 


vance, 

For all nonlinear elements falling into the sub- 
rroup defined in Section 6, the current i is deter- 
lined uniquely by the voltage v, and v in turn is de- 


sermined uniquely by the set of the V;, 9;, #; (1° 


Se 2s en). G; and C; can then be expressed formally 
1S 


G; = Gj (Vj 594,43) Nag eB CT 34) 


C; - C(V5 94543) Tee hall, (7.5) 


‘he G, and C; are continuous functions of this set of 


B iehles Eeaiutete except possibly at points V; mi) 
here they may have a pole of order one. 
_ Within limits of observational error, it is always 
-ossible to find equations of the form 


Py] tTPoWote** tphYy al (7.6) 


here pj,P9,°**p, are integers. The values assumed by 


5 and C;; as the sum ZI px | becomes very large are 
fined as "basic equivalent conductances and capaci- 


ances," and represented by the symbols G! and C'. 
J 
nerally the prime is used for all symbols associated 


ath G' and C'. It can be shown that it is then pos- 
; J 2 
ible to adjust the phases 91,99°**9, close to any 


lue by a simple shift of the time axis. On the 
ther hand it is quite clear from physical considera- 
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tions that such a shift in the time axis, which is a 

purely mathematical operation, cannot influence the 

values of G' and C'. This means that the basic equiv- 
J 


alent conductances and capacitances are independent 
of the phases P11 PQ°°* Op- G' and C' can then be ex- 
J J 
pressed formally as 
2 = iietkers, Aol ocean Cha 


Cc! 
J 


C1(V;w; ) ii Sie2ae bene (7.8) 
J 


The symbols G" and C" stand for the expressions 
J} 


G" = G; ao G! (7.9) 
J d 

Ce. 2 C; ey A (7.10) 
J n) 


and generally the double-prime symbols are used for 
all quantities associated with G" and C". If all 
J 


J 
=| p;, | become very large, then G" and C" approach zero. 


Equations (7,9) and (7.10) can also be written in 
the form 


G; pads 4 (7.19 
J i 


C; Sr! ot GS (7.129 
J J 


Then G' and C' express the properties of the non- 


linear element.if there is no synchronization or sub- 
harmonic resonance in the circuit of which the non- 
linear element forms a part. They can be measured — 
with relatively simple equipment; or if the internal 
mechanism of the nonlinear element is known, it is not 
too difficult to calculate them. The G" and C" have 


to be calculated separately for each specific set 0: 
P1»P9,***P, and cannot be measured easily. 
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Wave Wave 
Ana/yser Analyser 


Figure 7.2. Circuit for the Direct Measurement of G! 
je 


In many practical cases the C' can be assumed to be 
J 
zero. It is then particularly easy to measure the G!, 


J 
Nith a circuit of the form shown in Fig. 7.2. In 
this circuit the voltage v is generated by oscillators 
in series with the nonlinear element. The amplitudes 
»f the components of current and voltage, I; and Vis 


at the angular frequencies W; are measured with a 


wave analyzer or a similar instrument. By definition 


she value of the equivalent conductance is GinLi/V;. 
-t_can be assumed that all = |p,,| are large if none of 


che oscillators are synchronized. 


ANY 


oe 


Q / 2 3 
Y% 


Figure 7.3. Diagram of Constant Gi 


as a Function of Vy and V5 
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For a circuit having only two components of volt- 
age, the Gt can be plotted in the form of constant 


contour lines as shown in Fig. 7.3. 


8. Static Current-Voltage Characteristics 


In some special cases it is possible to calculate 
the G; and Cj. In the most important of these the de- 


pendence of the current on voltage is given in form of 
a static characteristic which is independent of fre- 
quency and has a fixed operating point. This type of 
dependence is shown in Fig. 8.1. It 1S important to 
realize, however, that such a characteristic is always 
based on a simplification and that in most physical 
cases the position of the operating point on the char- 
acteristic will change if the amplitude of oscilla- 
tion is changed. 


Figure 8.1. Static Current~Voltage Characteristics 


With this type of characteristic, then, the instan- 
taneous current is a unique, continuous function of 
the instantaneous voltage v; that is, 


i, rye Too a) master) (8.1) 


ce) 


where x is defined as v/V,, and I, and V, are con- 


stants with the dimensions of current and voltage. In 
particular V, is chosen so that -1$xs+l. This is done 


to permit the use of some functions which are defined 
only in that interval, as for example Legendre poly- 
nomials, It can then be shown that this nonlinear 
element falls into the subgroup defined in Section 6 
of this chapter. 


‘ 
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By finite sums of analytic functions, such as power 
Series or orthogonal polynomials, 9(x) can be approx- 
imated uniformly in the domain -l¢xs+l, or 


m 
O(x)¥ 2 ayyy(x). (8.2) 
‘el 


As indicated in the introduction to this chapter, 
the voltage v can, in all practical cases, be approx- 
imated by a finite sum of trigonometric functions. 
Consequently x is 


n 
> hs Re X;cos(wjtto;) (8.3) 
b ak | 
where 
\. 
X; =_ J 
V6 
also 
n 
ple It 2 I, jcos(wjtto;) tI, jsin(wjtto;) 


eee! 
+ currents at combination frequencies 
and higher harmonics. (672) 
The components of the current i at frequencies W),Wo, 


***W, can be calculated by the integrals 


— 
{ 


dee tinates poik, Hee: 
ej Te FE ier 1 cos(w itt; )dt (8.4) 


: + 
Tu = lim 1 fF ie sin(W;tto;)dt. (8.5) 


c ‘sf, 
UIT ey Sada 


Equations (8.4) and (8.5) can also be expressed in 
operational form as 
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Igj * {Py,i) (8.6) 


Ti -- waa 344) (8.7) 
J 
It can be seen that these operators are linear. 


Using them the equivalent conductance and capacitance 
can now be expressed as 


Gj ° TV, (Py id= 7 7, Pugniey? 
G, 
; x; { Py, ola } (8.8) 
C, -- —— Cups SOUARY Leakiear 
sejoued Mies Mabeat if 
_~— POLerrip Ca (8.9) 
i J 


where G=I,/V, is a constant with the dimensions of a 
conductance. 


9. Series Expressions for the Equivalent Impedances 


As shown in the preceding section, g(x) can be ap- 
proximated by the series 


m 
g(x) ¥ 2 ayy (x). (8.2) 
NZL 


Introducing the right-hand side of this equation into 
Equations (8.8) and (8.9) and making use of the lin- 
earity of these operators gives 


a, ta } 
G: = —{P,.°9(x)}x—— 2 aiP,.*v,(*) 
I X; si | X; AZ=1 US ee uN 
m 
vy 2 aGjp (9.1) 


AMI 
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m 
G;¥ 5 ayCiy (9.2) 


vhere 


Cet eer 1 Qe WR UX) Di 
jh ~ “Xj; T2057 VA 


For a given set of yw), the first step towards the 


-alculation of the equivalent impedance is the calcu- 
.ation of the Gi and Cin: For a specific voltage- 


‘urrent characteristic, that is for a specific set of 
i,, the equivalent conductance and capacitance, G; and 


4 can then be determined rapidly by adding the dif- 


ferent components. [or such a calculation, it is nec- 
>ssary that the a, are uniquely defined and become 


-apidly smaller with increasing A. 


Figure 10.1 Current-Voltage Characteristic with Spike 


0. Operational Calculus 


The function 9(x) has to be approximated by a fin- 
ite sum of analytic functions in such a way that the 
yerformance calculated from this approximation ap- 
proaches the performance based on (x) itself. The 
approximation in the mean is not necessarily the best 
way of doing this. For example, in the characteristic 
shown in Fig. 10.1, the spike in many cases will de- 
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termine the behavior of a system containing such an 
element. An approximation in the mean will, however, 
not adequately express this spike unless a great num- 
ber of terms are considered. 

It is important that the series converge rapidly 
and that the coefficients be uniquely determined and 
independent of the number of terms considered. Such 
series can be obtained for small neighborhoods of the 
origin by a Taylor series expansion and for larger in- 
tervals by orthogonal polynomials. These two cases 
are treated in some detail later. It is felt, how- 
ever, that these two methods are not necessarily the 
best and that other ways of approximating the charac- 
teristic (for instance piecewise approximation by 
simple analytic functions) might, in many cases, re- 
sult in more concise expressions. 

For both power series and orthogonal polynomials 
the basic terms are 


{P, *x\} and (Q, 0x1) (10.1) 
J 


In this section a method is discussed that will permit 
a, rapid calculation of these terms. 
As a first step an additional operator P, is intro- 


duced that determines the "mean value" of a time-vary- 
ing function, or 


lim 1. 
{P,°9} = Tos oT ie T edt. (10.2) 


Because of its simple structure, P, is particularly 
well suited for calculations. It can be seen that 


Drees, nice At+1 
{P Pos” } Mi OX; {Pi cigeht ht (10.3) 
and 
ae iia At1 
{Qu xt} in OFX, 5 caLD nel She (10.4) 


With the aid of these faeany the calculation of 
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P,.*x)} and {Q,,.°x"} can be reduced simply to that of 
J J 
Porxd* hy, 
The case where x has only two components 
> X} cos()t*9]) + Xgc0s(WottPo) (10.5) 


s discussed in some detail because it is of parti- 
ular importance. The frequencies W, and Wo are re- 
ated by the equation 


PY) ~ qWo. (10.6) 


he higher powers of x can be expanded in a binomial 
eries 


At] 
{P,-x\*1} = X; “{P,-cos**1(w t+p1) 4°: 
At] At+1L-k Atl-k 
+( re )Xy XM{P,*cos (w, t+9]) cosk( wottoo)} 


al 
#2294Xq  {Po*cos**](uyttyo)}. (10.7) 


he calculation of {P,+xh*1} is thus reduced to that 
f the terms 


: {P,*cos**1-k( uy t+91) cosk(wot+yo)} Kke0 S1aceet bs 
(10.8) 
The next step is to expand all terms of the form 
os\*1-k(y, t+9,) and cosk (Wott) into Fourier series. 


he coefficients Be ar ioe By of the Fourier series 
cos’a@ = a+ 2 a: cos ja (10.9) 


1) 
jr? 


ave the values 


v 
ay sere Hq if v>0 and even (10.10) 
2 
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=) if v<0 or odd 


= 1 sc 


if (v-j)20 and even (10.11) 


Si} if (v-j)<0 or odd. 
Using these relations, (10.8) can be expressed as 
{P* [b,+2b; cos j(wjtto7)] le, *2e; cos j(WgttPo)} 
(10. 12) 


The only terms contributing to the mean value 
(operator P.) are those independent of time. One such 


term is clearly b,c,, but a number of additional terms 
can be obtained because of the relation pw )*qwo. In 
fact all products 


byp cos ap(w tt, ) Cyq COs aq(WottPo) 


b 
= Papfaq {cosa( pp }-q9o) tcos(apwy ttaquottapp;taqpa) } 
2 
(10. 13) 
yield a term 
ba Sag 
re (App, -AqP9) (10. 14) 


which is independent of time. With the use of these’ 
relations, the expression (10.12) can be written as 


ee | 
Poo"? * Pap aq cos (4p 1-Aq9o)-. (10. Tae 


The terms b,c, and bupfag differ from zero only if 


the following conditions obtained from (10.10) and 
(10.11) are satisfied 


bocoF 0 if A+l and k are even (10,16) 
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baptaq F 9 if (A+1-k-ap) 7 0 and even 


and (k-dq) 2 0 and even.(10. 17) 

Adding the last two inequalities, the necessary con- 
dition 

Chee | - Of = dq) < 0 (10.18) 


is obtained. From this it follows that 


bapcaq = 0 if A+ < a(ptq) © (ptq). (10. 19) 
Therefore all products bapfag are zero if (ptq) be- 
comes very great and in particular if it becomes 
Msreater than (A*+1). 


As stated in Section 7, the prime indicates that 
(p+q) is very large. If the series is broken off 


after the term era and (ptq)>m, then all terms bapfag 


are zero. The only terms remaining then are the b,c, 


Pee cos M*1-K( 5, t+) cosk(wottog)} = boe 


oO 


k 
- " At+1-k Ae (10.20) 
2 


I 
/ 


The operator Pj’ is defined as 

(Pls) = {Py} - (Pte) (10. 21) 
therefore 
ae cos Mek Ca) Dy ) cos K(wottQo)} 
i byp lag cos(ap),-aqo). (10. 22) 
As an example, the term 


{pls xP} 
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is calculated. From inequality (10.19) it can be seen 
that all terms with exponents lower than (ptq) are 
Zero 


{Pilex PEASE} sea 0 nines 102 eae (10. 23) 


Hence (A+1)=(ptq) is the lowest exponent for which 
this term differs from zero. The number @ for which 
inequality (10.19) is satisfied is a=]. Now that p, 
q, ® and (At+1) are known, the values of k satisfying 
inequalities (10.16) and (10.17) have to be found. 
For a=] and (A+1)=(ptq) these inequalities are 


(q-k) 2 0 and even (10.24) 
(k-q) 2 0 and even. C1025) 
The only value of k satisfying these inequalities is 
kG, 97 bus 20) 15. Lound tune 
f 
{P,*xP*4}=(P XT XG{P,* cosP(w tp) cosT(wytty)} 


_L ptq 
ae )X7Xgbp Cqgcos( ph -d%q) 


q 
ae (P79) XX, cos( po} -a99) 


1 
=ohC MT) XBXJcos( p94 -499) (10. 26) 


Although this term is of considerable importance 
for phenomena connected with synchronization, even 
more important for the present discussion is the 
calculation of 


{Plena sb} = XM UPS cosh (wy t+91)} Te wit 


Nee Ab bek i 
oe k abe Xo Awe os\*1-k(w, +1) cos*(wottoo)} 


A+ 
°° 4X5 “{P!*cos**l (wottpg) }- (10. 27) 
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The general term of this binomial series is 


Pala thle kes i At+1-k 
( ‘i )Xy Xo TP rCos (w) t+o1 )cos*(wot*o5)} 


APLOS 1 kok 
Bee eh) a boc 


Oo-O 


A+1 i; X} Xo if (A+1l) and k are even 


gt £33 
2 512 (10. 28) 
=) -+f- eitherst Atl); on k sais: odd. (10.29) 
Hence 
‘ AFL 12 
ae ad CA Bane SE AT POS eeAFI-K Shs 1 
PS OEE eae SRL B.S Poa 70d BOR ts ao ee 
- AS 
if A’ is odd (10.30) 
= 0 if \ is even. (10.31) 


The usefulness of this result becomes evident in the 
next section. 
ms POWER -SERIES 

The current i=I,9(x) can be approximated by the 


rapidly convergent power series 
| m 
i¥T, 2 ayx' -1sx$1 (11.1) 
1 Set 


in certain cases where, for example, the expansion is 
imited to a small neighborhood of the origin or where 
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the voltage-current characteristic deviates only 


slightly from a straight line. The main advantage 


of 


this representation by a power series is the ease with 


which the equivalent impedances can be calculated. 
With vy=x" the considerations of Section 10 can 


applied directly 


G 
ea Mm «26 3 A+] 
G5. = tBu hh) eee ee ab 
J A as (+1) X;O%5 % 
G 
re) 2G, ) 
C-, = - {Q. ch re ee 
yA Pe Palatal ie Die ba eg 
Xj; j (A+1)X; oF 39; of 


The equivalent conductances and capacitances can 


obtained by adding the G5) and Cir 


m G, m x 
Gs =" 2 fatGay “hee Pe’ 2 cays 
MMi = Xe nN 
4] rN21 J 4 a Azl1 
2Ge fe m ay xh*1 
2 ' -{P ; 2 At] } Glas 
X; OX; re) oa 
m fe m i 
C) > a, C-y = 2{Q he aX } 
Az1 J X a i A=] 
re) m REL 
ay ee 2G, [pans a is Clin 
X.2y. Oj ° r=1 AFL 
From Equation (11.4) it can be seen directly that 
aG: AG: 2 At] 
Xj aol Gat 0G ee 
X; X; Xj X; . AML. okST . CEL; 


This result is used in Chapter IV. 


be 


2) 


3) 
be 


4) 


5) 


6) 
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Another important relation can be obtained from 
Equation (11,5) 


n n m 
2 X;20j2C; = - GE X;0{Q,.* 2 a,x} 
ay ae J r=1 


m 
+e**+X ww sin(w,t+9,) | a a,x dt 
A=1 
+7 A+ 
- G 1 i ees . clas 8 
OT. T +T dt A=] A+1 


bes ay x 


Aa +7 
cnligsee 
x Rea =] A+] age 


aah, bis ay 


In particular, for a system with only two degrees 
»f freedom 


X17012C) + Xy2uW2Cy = 0. (11.8) 


[t 1s thus only necessary to calculate one of the two 
>quivalent capacitances, the other can then be ob- 
tained from Equation (11.8). For two degrees of 
freedom, the values of Gin and G5 are given by Equa- 


cions (10.30) and (10.31); they are 


At] X+ At] 2 et 
(AF1)X; 2 9 J 1 
At] 
sath der Xo ] for \ odd 


(11.9) 


at for \ even. 
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For specific values of \ the Gin are 


Gip~0 Go9°0 
Gi 1°G, Gp 1°, 
Gj 970 G5970 
Gj 3°26, (xq 2#2X52) G5 3°SG, (Xq?#2X)?) 
Gj 470 240 


Gy 5 2G, (Xy #*6X1 2Xp7+3X94) Gy 5°2G, (Xg#+6Xq2X)2+3X,4) 


Of some importance are the terms, Gin and Ch): where 


(p+q)=(A+1). They can be obtained directly from 
Equation (10.26). 


2G ) A+ 
‘ee — is) ier l 
1A Tax, ox fea! 


ios) AnL pt xPo2 
gr-1 q )XF] 


Xo Icos( po) -q%o) 


G 
ee Ay PZ 
sgh=1'q)*4 X94cos( po] -qPo) (11.10) 


and similarly 


p-2 : 
Cis -4.0. (hixeiaxgdean Goat cee) (Listh 
lA gh-1, q 4 24 1 2 


From Equation (10.23) it follows that 
Gia Gri 0 for \<pt+q-l. (11.12) 


For each separate set of p and q the Gh and Ci have to 
be calculated separately from Equation (10,22). 
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In the following chapters the current-voltage 
characteristic is often given in the form 


a5 avtBv2trvstév4tevotess, ee ee 


where it 1s desired to obtain the results in terms of 
the amplitudes V|,Vo9,°°*V,, and the parameters a,6,°*° 


e,*°*’*. This can be done with the use of the following 
relations 


ea 
aR 
_ £04 
a = Gay Sa aac 
(@) 
ole G,a9 
V 
(@) at ee 
N = G,a3 be 
Vo2 


In the following example the equivalent impedance 
is calculated for a current-voltage relation 


i = 1,(x-0. 1x2+0,01x3-0.001x4) 


and a relation 2W 1 =3W9 between the frequencies (p=2, 
a=3). The Gy}; Go then have the values 


G} = Go[1+0.0075(X,2+2X52)] 


QQ 
No- 
Ut 


G,(1+0.0075(X_*+2X17)]. 

yf all the p's xAt1 only pie x9 is not zero because the 
condition 

: NeAeL cep tq 


_s satisfied only for \ = 4. Hence 
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G 
0.002 Yo O 

Nos 2 rip ade eae 8 1, 5 
G! = -0.001G%4 meme 

J J 
G 
0.002 “o 6) 

tae th = : datah Loess 

Ch = -0.001C}4 = -- x2 B95 {pie. 


The only value of k which satisfies the conditions 
At+|-k-p 2 and even 
k-q 2 and even 


is k=3, and the only product bap fag that is different 


from zero 18s agb3, therefore 
{Pls x9} = (3 )X1 2X9? {P,* cos*(wy t+01) cos?( aot top) } 


| Aap : i 
a 5(3)X] 2X52 8 cos( 291-399) 


5 
= 9X1 °Xq%cos( 201-399) « 
The values of Gj, G5, and C} are obtained by differ- 


entiation with respect to Xys X9 and ?] 


25 
= oe G,2X93cos (201-305) 


= -0.0005 GX5%cos( 291-399) 
0.002 


= -0.00075 oe 


0) 0022 Saran 3 

Cc! = G, 22 sin( 201-395) 

1 17392 
Aste 8 Wy 
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X,2 
= 0.0005 G,—% sin(29}-305). 
va 
These values could also have been obtained directly 


from Equatyone Ct bid)" and “CLT. P1)% ) can be ob- 


tained from the relation 


cy = -<} 


XOX 


u 


-0.00075 G, sin(291-399). 


a 


The total expressions for the linearized impedances 
are, therefore 


= G [1+0.0075(X12+2X5") -0. 0005X%cos( 291-399) ] 


op) 
— 
I 


Gy = G,[1+0.0075(X_2+2X17)-0.00075X)7X_cos( 291-399) ] 


X.° 


: C) @10) 005-6... sin(29)-399) 


Wy] 
XaeX 


Cy = -0.00075 G, sin( 291-309). 


“ 
As was mentioned before, the current 1=1,9(x) can 


e€ approximated in a domain -1sx$l by a finite series 
of orthogonal polynomials. The coeffieients of this 
series can be obtained from orthogonality relations. 
As a rule they will become small rapidly, even if the 
curve differs considerably from a straight line. This 
inite series of orthogonal polynomials can be trans- 
formed into a power series, but the coefficients will 
not necessarily converge and they will also depend on 
the number of terms of the orthogonal lines. 

The simplest way to calculate the G5 is to break 
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up each individual polynomial into a sum of powers of 
X, to operate separately on each of these, and then 
put the individual components together again. 

As an example, the Gi are given for Tschebyscheff 


polynomials and a voltage with two components. 
Hera N= 2,4,6°°° 
Giy ~ G, 
Giz = 3G,(X12+2X52-1) 
G}g = 5Go( 2X) 4+12X12X52+6X54-3X) 2-6X52+1) 
Giq = 7G,( 5X, °+60Xy 4X52+90X) 2X54+20X5%10X, 4 
- 60X] 2Xq%-30X94+6Xy2+12X52-1) 
The values for G5 can be found by interchanging the 
X; and X59. The ratios Gj} 3/G, Gis/Go, and G}7/G, have 


oY & 
ee 
8 


tees 


%, Ariplitude of Oscillation at Freguenc 


0 02 0. 
x,, Amplitude of Oscillation at Frequency a, 


Figure 11.1. Curves of Constant Gi3/S as a Function of X-Xy 


%, Arnplitude of Oscillation at Freguercy Wz 


x, Anolitude of Oscilatiora at Freguericy Wz 


BUL. 395. NONLINEAR SYSTEMS 


0 02 08 \\ 
x,, Anyplitude es peepee at Frequency %& 


Figure 11.2 


Curves of Constant G*_/G. as a Function of X,—X 


SVEN; Lee 


= 


TM A | 
JAS) 


of 
A 


_———— 
ma et 
_ 
ay 
ey 
0 
0 02 OF 06 08 1.0 
%,, Anmplitude of Oscillation at Frequency w, 
Figure 11.3 


Curves of Constant 47/6, as a Function of X,-X> 


oe 
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been plotted in Figs. 1101," 11,2; and “liv one ait can 
be seen that these values are relatively small for 
large portions of the X)-X» plane. 


The calculation of the G; has to be carried out 
in two steps. First, it is necessary to determine how 
many terms of the series need to be considered. This 
can be done, for example, with the help of Figs. 11.1, 
11.2, and 11.3. In a second step, the sum of poly- 
nomials has to be converted to a sum of powers of X 
and from here on the procedure is identical with that 
used in Section 10. 


lll. HIGHER APPROXIMATIONS 


The discussion in this chapter concerns a method 
for obtaining higher approximations to the solutions 
9f non-linear differential equations. The method 
nakes extensive use of the physical concepts under- 
lying the problem. 


| 2. Introduction 


The differential equations occurring in nonlinear 
nechanics can not as a rule be solved explicitly, and 
it 1S consequently necessary to approximate their 
solutions by sets of analytic functions. This approx- 
imation of the solution by a set of analytic functions, 
»referably orthogonal, is often used in mathematics 
for the solution of such problems as linear differen- 
tial equations. [For them it can be proved that under 
certain conditions the set of analytic functions con- 
verge to the solution$’?’ No such statement, however, 
can be made regarding the periodic solutions of the 
original nonlinear differential equations; we know 
1ext to nothing about the convergence of successive 
approximations. For the mathematician this situation 
is highly unsatisfactory but for the engineer it 1S 
somewhat less so. In contrast to the mathematician, 
the electrical engineer is less concerned with the 
rigorous mathematical soundness of his calculations as 
long as they give results that confirm experimental 
neasurements and allow a certain amount of extrapola- 
tion. 

In many physical nonlinear systems (except, for 
=>xample, those having highly saturated iron cores), 
the first approximation gives results which are AS 
quate to explain the phenomena qualitatively and which 
also correspond well to experimental measurements. 
[Therefore a first order approximation is entirely 
satisfactory in most cases. In some instances, how- 
ever, where marked discrepancies with the experimental 
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results occur it 1s necessary to consider also the 
second and possibly even higher approximations. 

For autonomous systems having only one degree of 
freedom and, therefore, permitting a representation 
in the phase plane, it is possible to check the qual- 
ity of an approximation by comparing the phase trajec- 
tories corresponding to the analytic solution with 
those obtained by the topological method of isoclines 
Ciigs el? Ayers) = 

In these figures, the limit cycles corresponding to 
periodic solutions of the van der Pol equation, 


ee 


U - p(l-u2)a + u = 0 eee 


Figure 12.1 
Limit Cycles of the Van Der Pol Equation for U-0.5 
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Figure 12.2 
Limit Cycles of the Van Der Pol Equation for U=1.0 


ave been plotted in the phase plane (abscissa u, 
dinate w=u) by the method of osiclines for u=0.5 
ad w=1. Also in the same figures the function 


wei 2ircds* t -— sin 3t C122) 


as been plotted. As is shown later, this function 
»rresponds to a second-order approximation to the 
riodic solution of Equation (12.1). It can be seen 


iat for values of Lup to 0.5 this approximation is 
41te good. 
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The method used in the following is essentially 
that of Kryloff and Bogoliuboff. The main change is 
that physical reasoning has been substituted for math- 
ematical considerations. This change has helped to 
reduce the number of calculations by a considerable 
amount. 


13. Outline of the Method 


As a first step towards the approximate solution of 
a nonlinear differential equation, an electrical cir- 
cuit corresponding to this differential equation needs 
to be found. For a mathematician this step would only 
complicate the problem but for an electrical engineer 
trained to think in terms of voltages, currents and 
impedances, it puts the problem on a more familiar 
basis. Once this circuit has been found, the reason- 
ing can be’confined to the physical concepts of the 
problem. 

Instead of giving a general discussion of the 
method, a simple example is treated which illustrates 
all the important points. Problem: [Find the frequency 
of the periodic solutions of the differential equation 

U + wo(atbutcu2)a + w2u = 0. Bee 
It is known that this frequency 1s approximately 4; 
in the following development a second approximation is 
found. 

It is assumed that the solution of Equation (13.1) 
is stable; this means that a<0 and c>0.‘23) For b=0, 
this equation can, by a linear transformation in u and 
t, be transformed to the van der Pol equation 


u - u(l-u4)a + u= 0. (1.1) 

A circuit corresponding to Equation (13.1) is 
shown in Fig. 13.1. In this circuit the inductance L, 
and the capacitance C, form ‘an oscillatory circuit; NL 


is the nonlinear element; the instantaneous current 1 
is a function of the instantaneous voltage v given by 
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i = w(v) = av+Bv2t+yv3, (1329) 


11S particular function has been chosen in order to 
ytain correspondence with Equation (13.1). 


Figure 13.1. Circuit Corresponding 
to the Differential Equation (13.1) 


The differential equation of the circuit shown in 


g. 13.1 is 


Betsey a ae a 13.3 
Odt (avtBvetyv?) L Gy toe 0 ext 


oO 


-, by differentiating with respect to time and 
viding by C: 
es ; l 
yt — (at2RBv+3yv2) ¥ erste sO Gb3 hd) 
Ca LC 


oO 


uation (13,4) corresponds to Equation (13.1) if 


i 
ieee 
W) 
e LC, 
a = W Coa 
ee 6° Cb (13.5) 
2 WoC, ih 
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Thus, it is seen that the circuit of Figs 913, 1) does 
correspond to Equation (13.1) and further treatment 
of the example is confined to this circuit. 
First the voltage is assumed to be purely sinusoi- 
dales that yas 
Vo TAY) Cosmm es (13.6) 
This assumption can be justified by physical consider- 
ations if the current passing through the nonlinear 
element is small as compared to that in the inductance 


and capacitance. Then the instantaneous current 1, 


corresponding to this voltage can be calculated from 
Equation: (13:2). 


Second the voltage v, generated by the current i, 


passing through the linear circuit formed by L, and C, 
must be found. The voltage v; will be somewhat dif- 


ferent from v,; in particular it will contain higher 


oO? 
harmonics. Again a current i; 1s calculated from 
Equation (13,2). The amplitudes of oscillation and 
the frequency can then be found by replacing the non- 
linear element with a linear impedance according to 
the theory of equivalent linearization. 


14. The First Approximation 
Using the voltage 


vo = Vocosut (13. 6a 
the current es from Equation (13,2) is : 
pe avotBv,?tyv,° | 


8 3 
“QVo°t(aVo*71Vo%) cosut 


ib 
AV, 2cos2utt WV, Scosdut. (14.0 
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Figure 14.1. The Equivalent Linear Circuit 


It has been shown in Section 2 that to a good ap- 
roximation the nonlinear element can be replaced by a 
inear element that will produce the same components 
f the current 1 as the nonlinear element, but no 
igher harmonics. If the components of 1 at frequency 

are 


1/9 I) coswttlosinwt 


+ d-c term and higher harmonics, (14,2) 


hen the circuit of Fig. 13,1 can be replaced by the 
mnear circuit shown in Fig. 14,1. In this figure the 
onlinear element has been replaced by a conductance G 
n parallel with a capacitance C. The values of the 
equivalent" conductance and capacitance are 


G = nt 
ie (14.3) 
1 
D 
Cc =—., 
WV, (14,4) 


As can be verified easily, the current flowing 
hrough this equivalent linear network for a voltage 
> Yocosut is 


OG 1, coswttIosinut. MS EACSS) 
nis current has the same components at frequency W as 


he current flowing through the original nonlinear 
ircuit (Equation 14,2). 
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Comparing Equations (14,1) and (14.2), it is found 
that 
3 


: Ee 
es (14.6) 


ee (14.7) 


consequently substituting these values for I, and I» 


into Equations (14.3) and (14.4) gives for the equiv- 
alent conductance and capacitance 


3 
Gita bapatilles (14.8) 
C280. (14.9) 
The frequency of the equivalent circuit is 
Pole op 
a) SLLEGN a Wo (14. 10) 


For steady-state oscillations, it is necessary that 
G20, vores that 


a 
ele ac eee ey (14.11) 
3 Y 


which for the van der Pol equation reduces to the 
well-known result that 


V eres (14, 12) 


15. The Second Approximation 
Using the current 1, from Equation (14.1), the 


voltage vj generated across the linear L,C, network is 
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- —————— sin3ut. ULSS)) 


In the first approximation ® is equal to w). 


therefore the voltage v); is 


Ae V,coswyt- (EV,2sin2ugt+ WV,3sin3ugt) 


WC 


o~o 
= b ed iy 3. 
= NSE rd fe SAAS teers fe sin3u,t Ch5.2) 
yhich for the van der Pol equation reduces to 
= Qcos t - — sin3t (12. 2) 
i Giada -—sin3t. : 
l 4 


‘or a given V, the amplitudes of the higher harmonics 
ire small if the product w,C, is large. [or a phy- 


sical oscillator, the damping of the oscillatory cir- 
sult 1S not zero, and it can then be shown that the 
implitudes of the higher harmonics are inversely pro-, 
rortional to the quality factor Q of the circuit. 

The current 1, can now be obtained by using Equa- 


lon (13,3) again. Since the approximation is not 
-arried out further, it is not necessary to calculate 
‘very term of 1). The correction in frequency is de- 


ermined by the equivalent capacitance C and C in turn 
s determined only by Ig (Equation 14,4). Conse- 


uently only the terms of 17 contributing to I5 are 
eeded. They are 


0340 V,2sinQut coswWt 


suc Vo™ cos2ut sin3wt CLS5.3) 
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from which the value of 1) is 


l 2 ie 
_B* he Pe v5, (15. 4); 


I 2 
2 l 9 
(ete ee B V,2 Praag of (15364 
WV, 3 w2C, 128 w2C 
dl eS Py p )2 oe 4 9 i )2y 4 
C, BOG 128 wC,, ° 
t l 
aise 
= vV,2 + ; c2V,4. (15.6) 


The oscillatory frequency of the equivalent cir- 
cult shown in Fig. 14,1 is 


Dea nié ged 22a san lien 
Ovi CatG) ano aC) ee 


It2 


Le iA 
AS = 1b, Be ZV. 4, 15% 
Wo (1 2 C, erik! oA appv: yo eal ee 


From Equation (15.7) it can be seen that the oscilla- 
tory frequency is always decreased by the nonlinear- 
ity of the circuit. 

The expression for the oscillatory frequency ob- 
tained in Equation (15.7) could also have been ob- 
tained directly from Equation (69.16) in Minorsky. ed 
As pointed out at the beginning of the chapter, it 1s 
the method rather than the result that is new. 

The oscillatory frequency corresponding to the 
van der Pol equation is 

4 2 
Beek Tae (15.8) 
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the correction in frequency will be 


lor vacuum tube oscillators, Equation 
does not 


most cases, 


[uite small. 
15.7) holds-only if the operating point 


thange as the amplitude increases. 


IV. SIMULTANEOUS OSCILLATIONS 


This chapter concerns simultaneous oscillations ir 
nonlinear systems with several degrees of freedom. As 
shown by L. Skinner, 7° such simultaneous oscilla- 
tions are possible if the nonlinear element satisfies 
certain conditions. 


16. Introduction 


It is often assumed that an oscillatory circuit 
with several degrees of freedom, such as the one showr 
in Fig. 16.1 will not oscillate at more than one fre- 
quency at a time. As is shown in this chapter, suck 
an assumption is erroneous and it is possible for @ 
circuit with several degrees of freedom to oscillate 
simultaneously at several frequencies. This circuit 
consists of a nonlinear element NL and n oscillatory 
circuits formed by the inductances Loj) capacitances 


Figure 16.4 
Oscillatory Circuit with Several Degrees of Freedom 
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j? and conductances Go; - The instantaneous voltage 
‘ross the nonlinear element is v and the current 
issing through it is 1. 

A circuit similar to the one shown in Fig. 16.1, 
it’ with only two degrees of freedom and a nonlinear 
ement having a current-voltage relationship 


ie tke avtyv9 Glo.h) 


iS investigated by van der Pol in 1924. He found 
rrectly that such an element does not permit simul- 
neous oscillations in the circuit 

The calculation was extended in 1948 by L. Skinner 
. a current-voltage relationship “°) 


it avtBvetyyo+6vetev, CL652) 


‘inner proved that with this characteristic simul- 
ineous oscillations are possible. The calculation 
in be extended to three or more degrees of freedom 
id current-voltage characteristics of a degree higher 
ian the fifth, but then it rapidly loses simplicity. 
1 another paper W. Huggins investigated simultaneous 
scillations observed in reflex klystrons, 

The case where two circuits are synchronized with 
ich other and their frequencies have a fixed ratio, 
r example, 2:3 has been treated by R. Fontana, ‘7?’ 


'. Systems with Two Degrees of Freedom 


The object of this section is not so much to give 
general discussion of simultaneous oscillations as 
; 1s to show that such oscillations are possible. 
le existence of these oscillations are demonstrated 
‘rr the circuit with two degrees of freedom shown in 
g. 17.1. The two oscillatory circuits are formed 
* the capacitances Coy» Cog, the inductances Lp}. 


2 and the conductances Go1) Goo: It is assumed that 


ae quality factors of both circuits are large. Fur- 
‘ermore, it is assumed that the resonant frequencies 
’ the two circuits are not too close together. 
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Figure 17.1. Oscillatory Circuit with Two Degrees of Freedom 


It is then possible to approximate the instantan 
eous voltage v across the nonlinear element by 


Va V1 cos(#1t*91) *Vocos(WottGo) Olt. 
where V, cos(,t*+)) is the voltage between point 


A and B, and Vocos(Wott9o) between points B and ( 


The angular frequencies , and Wo are approximate] 


equal to the natural frequencies of the oscillator 
circuits. 

Within limits of experimental error, it is alway 
possible to find between the frequencies , and Wo 
relation 


PW, ~ qWo (10.6 


where p and q are integers. 

The case where p and q are small integers 1S simi 
lar to the synchronization of an oscillator by 
external voltage, with the difference that the two os 
cillatory circuits synchronize each other. A ver 
good treatment of synchronization has been given bh 
Fontana.'?*) The "synchronizing force" becomes ver 
small if the sum (ptq) is large, and the synchroniza 
tion terms can then be neglected. Skinner assumed 1 
his paper that W 1/9 was an irrational] number in whic 


case the sum (ptq) is then infinite. The same assump 
tion 1s made in this chapter. 
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.» The Linearization 


It 1s assumed that the voltage v across the non- 
near element and the current i passing through it, 
e related by Equation (16.2). Using this equation 
r the voltage v, this relation can be expanded in a 
uble Fourier series 


3 
3 Vy [at-n(Vy?#2V592) 
#32(V M464 2¥92+3V94)] cos(wyt*o1) 


Vp [at2A(Vq2#2V 2) #20 Vo4+6V2vq2+3V4) ] cos wyt toy) 


+ d-c term, combination frequencies, 
and higher harmonics. Gigs tu) 


cause the sum (ptq) is infinite, no combination fre- 
encies or higher harmonics are equal to W) OF Wo. 


As has been stated in Chapter I, the behavior of 
e circuit will remain nearly the same if the non- 
near element is replaced by an equivalent linear im- 
‘dance so that the components of the current at the 
quencies W, and WwW. remain unchanged. Therefore the , 


irrent passing through the linear impedance is 


3 
q Vy la+Sy(V,242V92) #2e( Vy 4+6V4 2¥52+3V94)] cos (uy tty) 


3 
g lariy(Vo2+2V 2) +e (Vo4+6V, 2V52+3V4 4)] cos(uyt tog). 
(18,2) 
nce the voltage across the linear element is 
Vie Vj cos(W1tt1) tVocos(Wot*Po), GEL) 


_can readily be seen that this equivalent linear im- 
dance must be the conductance G shown in Fig. 18,1. 
1s conductance is a function of the amplitudes Vi: 
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Vo and the frequency w. Its value at a frequenc 
Wy] is 


3 5 
Gy = aty(Vy?+2V9?) + e( Vy HOV, 2VQ7#3V04) (18.2) 
and at frequency Ww» 1s 
= qty Vo 2+90V, 2) el VottEV, 2V0243V,4 
Go = Pern: TaN eee +6V1 Vo *3V1] ne (18. 3) 


The direction of current and the sign of the coeffi- 
cients @,***e has been chosen so that a positive con- 
ductance corresponds to a dissipation of power. 


Figure 18.1. The Equivalent Conductance 


At frequency #, the impedance of the upper circuit 
in Fig. 18.1 is very large and that of the lower cir- 
cuit comparatively small so that it can be neglected. 
Consequently the circuit reduces to Fig. 18, 2ag 
Similarly, at frequency Wo it reduces to Fig. 18. 2b. 
ln the first approximation the behavior of the system 
represented by these two circuits 1s now identical 
with the behavior of the original circuit 4 


aQ-At Frequency ©), i b-Ar Frequency Ge 


Figure 18.2A. Equivalent Circuit at Frequency W, 


Figure 18.2B. Equivalent Circuit at Frequency W, 
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). Equilibrium and Stability 
If the amplitudes of oscillation are to remain 
ynstant, then the damping of both the circuits of 
gs. 18.2a and 18.2b has to be zero; that is, 
m | 5 
3 1*Gy “Gy 1 *04 V(Vy2+2V9?) +0( Vy A+6V) 2V92+3V 54) =0 
(19.1) 
3 5 
30.2*Go "Gy tat y(Vo?#2V) 7) +2 (VoH+6V 2V9?+3V14)=0. 
(19,2) 


1 the following section the symbols Vo, and Voo are 
sed for the values of V, and Vo that satisfy both 


lese equations. 

The system is in equilibrium if Equations (19,1) 
id (19,2) are satisfied. This equilibrium may or may 
»t be stable, however. An.equilibrium position is 
sable if after a small disturbance the system returns 
) 1ts original state.. After a small disturbance, the 
iplitudes V) and Vo will have values slightly dif- 


srent from Vol and Voo 
Vo = Vo2 + 6Vo5 6Vo << Vo: (19.4) » 


is assumed that the disturbance occurs at the time 
0. At this time the values of 6Vy and 5Vo are: 
HAV: bVo=AVo. In this section the differential 


uations for 6V) and 6V5 are found under the assump- 


on that no additional disturbance occurs. If both 
‘| and 5Vo approach zero as the time increases for 


1 values of AV, and AVo, then the system is stable. 


| The average energy stored in the circuit of Fig. 
2a 1s: 


1 
Ey = 5 Coi(Vort8V))?. (19.5) 
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Twice during each cycle this energy changes from a 
electrostatic to an electromagnetic form and back 
but its total amount remains approximately constant. 
The first time derivative is 


, = 


This must be equal to the difference in average power 
created and dissipated in the circuit 


: (Vo 1*8V})? 

grr prc rae A 5 A Ca 
Vous oya (19.7) 
Voot6Ve 


The minus sign in Equation (19,7) 1s necessary because 
a positive conductance will dissipate power and de- 
crease the amount of energy stored. The sum (G)1*G)) 
can be expanded in a double Taylor series around VoL 


Vo2 


(Vp 1+6V})2 0G, 3G 


1 l 
+G +—6V,}t+—6Vo) tet]. 
2 UG GDy., SV, lev, Dp ee 


Voo Voe 


E,=- 


(19.8) 
The first term of this expansion is zero (Equation 
19,1). All terms of higher order in 6V) and 6V5 will 
be neglected. The following equation is then obtained 


d(6Vj)_ OG, OG 


For the circuit of Fig. 18. 2b the equation is 


d(5Vy)__ 8Gy_ OG : 
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e behavior of 6V) and 6V5 with respect to time is 


en determined by this set of two, simultaneous, 
near differential equations whose solution is 


Le 
SV, = Ae LY + Be 2” (19.11) 


t 
8Vy = Ce LY + De'2 (19. 12) 


ere XK] and Ko are the roots of the determinant, 


Vo, Gy Ne iat ecaeep 
2G), Vy, 2G), Vo 
= 19, 13 
02  0Gp Mie) ve Retcatier, pi cco ueaa 
oe yt Pipa gem she 


d the coefficients A, B, C, D, are determined by the 
itial conditions. 


| It was stated at the beginning of this section 
t the system is stable if both 6V) and 6V5 approach 


ro as the time increases. This happens if both Xj 
| Xo in Equations (19.11) and (19,12) have negative . 


al parts. It can be easily verified that the real 
ts of %, and %» are negative if the following two 


aditions are satisfied 


SAE OST AES 


EET (19.14) 
C1 Oy Cog Vo 
AGn OGs + hOGw 8G 
1 OG, 1 Ge 
— a IY G1O*15) 


V, OV>  BVy OV} 


For the special case where the current-voltage 
racteristic can be expressed by a finite power 
cies, it can be shown (Equation 11.6) that 
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OG, 9G 
Vol ay. < Vo2 av,” (19,1 
Zz | 
Using this relation, inequality (19.15) can be si 
plified to 
O62 OGy .u Vee 
1. 202. sy Yo ieee yas: (19.17 
QV, BV Von Vo 


It can be readily seen that the inequalitie 
(19.14) and (19,17) are satisfied if 


OG 

be >) (19.18 
1 

OG5 

mE tS 19,19 

BV, 0 ( 9 


ees Losin oan 
OV) Vo Voo OVo 


(19, 20 


Using the values for G, and Go from Equations (18,2 
and (18.3) gives 


OG 
ne 3.49 2 2 
OG 
we age 349 2 2 
oGe 
Abe 15 2 D) 
ay Vog[3y*—8(Vo1"*Voo )] (19.2 


so that the inequalities (19.18), (19.19) and (19.2 


become 


OY io 2 2 
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eee Vie ee3y 2)] > 9 (19.2 
02S Y*5°(%92°73%p1 -25) 
3 5 9 Sens) 5 


15 
> Vo1Voo(3yt-5eVy97)?. (19. 26) 


In order to simplify calculations, two new vari- 
les 
5 e€ 9 5 e€ 9 
ae ee Z% «== 19, 27 
ita, 0! meen 02 (19. 27) 
e introduced. Using these new variables, inequality 


9,26) becomes 


3-8(y+z) t6y2+6z2t+Byz < 0 (19, 28) 


ich is satisfied for values of y and z falling in- 
de the ellipse of Fig. 19.1. Since all values of y 
dz falling inside this ellipse are positive, it can 
en be concluded from Equation (19,27) that either ¢€ 
-y is negative. 

_ Similarly inequalities (19.24) and (19.25) can 
expressed in terms of the new variables y and z. , 
fo cases have to be distinguished: y>0 and y<0. 
ir y?0 the inequalities are 


1 - 3(y+3z) > 0 (19.29) 
bew 3( Gyr zy 0% (19. 30) 


se inequalities are satisfied in region I of Fig. 
.1. For y<Q they are 


Let aty tz) 0 (19.31) 
1 -- 3(3y*z) < 0 G19332) 


are satisfied in region II of Fig. 19.1. 

It can thus be seen that all three conditions 
9.18); (19.19) and (19,20) are satisfied if y<0, 

) and if the values for y and z fall inside the 
lipse of Fig. 19,1. 
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Figure 19.1. The Region of Stability 


20. Conclusions 


It has been shown that "asynchronous" simultaneou 
oscillations in circuits with two degrees of freedc 
are possible. As has been shown elsewhere, they mi 
indeed be obtained experimentally. An additional 11 
vestigation would show, however, that the oscillatior 
are not self-starting. This complicates the experi 
ment but does not make it impossible. 

Synchronous simultaneous oscillations (w1/W9=4/E 


where p and q are small integers) on the other har 
can be obtained experimentally very easily, but the 
are more difficult to study analytically. 

To complete the proof, it is necessary to sha 
that by a suitable choice of Go] and Go9 it is alway 


possible to obtain at least one equilibrium poir 
falling inside the ellipse of Fig. 19.1. Equation 
(19.1) and (19.2) transformed into the y-z system al 
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spe 

-Gyy = & + [-2( y+22) +y2+6y24322] (20.1) 
OL eS 2 2 

“Gyo = o an = — [-2( 2y+z) +z“t+6yz+3y~) . (20, 2) 


any arbitrary set of y and z there corresponds a 


tof Go 1 and Goo satisfying Equations (20.1) and 
0.2). Therefore it is always possible to adjust Go 
d Go9 in such a way that one equilibrium point falls 
side the ellipse of Fig. 19.1. 
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